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Abstract 

The investigation into the fast parton energy loss in cold nuclear matter is crucial for a good 
understanding of the parton propagation in hot-dense medium. By means of four typical sets of 
nuclear parton distributions and three parametrizations of quark energy loss, the parameter values 
in quark energy loss expressions are determined from a leading order statistical analysis of the 
existing experimental data on nuclear Drell-Yan differential cross section ratio as a function of 
the quark momentum fraction. It is found that with independence on the nuclear modification of 
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parton distributions, the available experimental data from lower incident beam energy rule out the 



incident-parton momentum fraction quark energy loss. Whether the quark energy loss is linear or 
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quadratic with the path length is not discriminated. The global fit of all selected data gives the 
quark energy loss per unit path length a = 1.21 ±0.09 GeV/fm by using nuclear parton distribution 
functions determined only by means of the world data on nuclear structure function. Our result 
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does not support the theoretical prediction: the energy loss of an outgoing quark is three times 
larger than that of an incoming quark approaching the nuclear medium. It is desirable that the 
present work can provide useful reference for the Fermilab E906/SeaQuest experiment. 
Keywords: quark energy loss, sea quark distribution, Drell-Yan. 
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1 Introduction 

The parton energy loss in high energy collisions has attracted an increasing amount of 
attention from both the nuclear and particle physics communities for over two decades. 
There is a rich theoretical literature on in-medium parton energy loss extending back to 
Bjorken, who proposed the suppressed production of particles having large transverse mo- 
menta, known as jet-quenching, was the "smoking guns" of the Quark Gluon Plasma (QGP) 
formation in high energy nucleus-nucleus collisions' 1 '. The wealth of experimental data on 
jet-quenching from RHIC' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ] and LHC' 8 l reflect clearly the energy loss of fast partons 
while traversing this hot and dense medium. However, a detailed understanding of the 
parton energy loss in hot and dense medium requires the good investigation into the fast 
parton propagation in cold nuclear matter because there are common elements between the 
two mediums. 

Two sets of experimental data from the semi-inclusive deep inelastic scattering of lepton 
on nuclei and the Drell-Yan reaction' 9 ! in hadron-nucleus collisions can provide the essential 
information on the energy loss of fast partons owing to multiple scattering and gluon radia- 
tion while traversing this cold nuclear medium. The semi-inclusive deep inelastic scattering 
on nuclear targets is an ideal tool to study the energy loss of the outgoing quark in the cold 
nuclear medium. In our recent article' 10 ', the experimental data with quark hadronization 
occurring outside the nucleus from HERMES' 11 ] and EMC' 12 ! experiments are picked out by 
means of the short hadron formation time. A leading-order analysis is performed for the 
hadron multiplicity ratios as a function of the energy fraction on helium, neon, and copper 
nuclei relative to deuteron for the various identified hadrons. It is found that the theoretical 
results considering the nuclear modification of fragmentation functions due to the outgoing 
quark energy loss are in good agreement with the selected experimental data. The obtained 
energy loss per unit length is 0.38 ± 0.03 GeV/fm for an outgoing quark by a global fit. 

The hadron-induced Drell-Yan reaction on nuclei is an excellent process to investigate the 
incoming quark energy loss in cold nuclear matter because the produced lepton pair does 
not interact strongly with the partons in the nucleus. A series of experiments' 13 ^ have been 
performed at Fermilab and CERN which presented the Drell-Yan differential cross section 
distributions in order to test the theoretical model, know the momentum distributions of the 
projectile and target quarks, and explore the nuclear target dependence. Four experimental 
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collaborations have measured the Drell-Yan differential cross section ratio of two different 
nuclear targets bombarded by the same hadron at the same centre-of-mass energy in order 
to study the nuclear effects on Drell-Yan reaction. They are NA3^ and NAlO' 15 ' Collabora- 
tions from CERN, and E772t 16 l and E866t 17 ^ Collaborations from Fermilab. The advantage of 
using the Drell-Yan differential cross section ratio at the same energy is that the differential 
cross section ratio can reduce the dependence on the beam hadrons, and cancel the most un- 
certainties regarding the lepton pair production. Additionally, the differential cross section 
ratio can avoid the influence of the QCD next-to-leading order correction. Theoretically it 
has proved that the effect of next-to-leading order correction on the Drell-Yan differential 
cross section ratio as a function of the quark momentum fraction can be negligible for the 
800 GeV proton beam at Fermilab and lower energy beam' 18 !. 

In our previous articles' 19 ' 20 ', the energy loss effect on the Fermilab E866 nuclear Drell- 
Yan differential cross section ratio was investigated as a function of the quark momentum 
fraction of the beam proton at the hadron level in the framework of the Glauber model, and 
at parton level by using two typical kinds of quark energy loss parametrization, respectively. 
It was confirmed that the energy loss effect can suppress evidently the differential cross 
sections versus the quark momentum fraction. In the recent work' 21 ', the study on quark 
energy loss is extended to the E772 data without performing the global fit to E772 and E866 
data. It is found that the quark energy loss effect on nuclear Drell-Yan cross section ratio 
becomes greater with the increase of quark momentum fraction in the target nuclei. The 
global analysis of nuclear parton distribution functions including E772 data overestimates 
the nuclear modification in the sea quark distribution if the quark energy loss effect is 
neglected. It is noticeable that the E772 and E866 Collaborations used the same 800 GeV 
proton beam incident on various nuclei. The measured momentum fraction of the target 
parton is in the range 0.01 < x 2 < 0.271. 

The main goal of the present work is to extract the incoming quark energy loss in cold 
nuclear matter systematically from a global analysis of these experimental results on nuclear 
Drell-Yan differential cross section ratio from NA3 [14] and NA10 [15] Collaborations at CERN, 
and E772' 16 ' and E866' 17 ' Collaborations at Fermilab. The main improvements over our 
earlier work are twofold: on one hand, the error estimate for the incoming quark energy loss 
is presented, and on the other hand, by adding the NA3' 14 1 and NAlO' 15 ' data from CERN, 
the used experimental data can cover 140 GeV, 150 GeV, 286 GeV and 800 GeV incident 
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hadron beam with the target-quark momentum fraction from 0.01 to 0.45. The extended 
beam energy and kinematic coverage significantly increase the sensitivity to incident parton 
energy loss and nuclear modification in the sea quark distribution. It is hoped to provide 
a good understanding of the parton energy loss in cold nuclear matter from the available 
data, and to facilitate the theoretical research on the energy loss of an incoming quark and 
outgoing quark in nuclear matter. 

This article is organized as follows. A brief formalism for the differential cross section in 
nuclear Drell-Yan process is detailed in Sect. 2, followed by the data selection in Sect. 3. The 
obtained results are discussed in Sect. 4. Finally, the summary and concluding remarks are 
given in Sect. 5. 

2 Dilepton production differential cross section in nuclear targets 

At the leading order(LO) in perturbation theory, the lepton pair production differential 
cross section in hadron-nucleus collisions can be obtained from the convolution of differential 
partonic cross section qq — > l + l~ with the parton distribution functions in the incident 
hadron h and the target nucleus A. With neglecting the incoming quark energy loss in cold 
nuclear matter, the differential cross section is written as 



where £1(2:2) is the momentum fraction of the partons in the beam hadron(target), a em is 
the fine structure constant, y/s is the center of mass energy of the hadronic collision, e/ is the 



variable x and photon virtuality Q 2 in the hadron (nucleon in the nucleus A), and the sum 
is carried out over the light flavor. 

In the hadron-induced Drell-Yan reaction on nuclei, the incoming quark can lose its 
energy AE q , owing to multiple scattering on the surrounding nucleon and gluon radiation 
while propagating through the nucleus. The energy loss of an incoming quark results in an 
average change in its momentum fraction prior to the collision, Ax\ = AE q /Eh, where Eh is 
the incident hadron energy. On the basis of theoretical research, three parametrizations for 
quark energy loss have been proposed separately by Brodsky and Hoyer' 22 !, Baier et al.' 23 ', 




(1) 



charge of the quark with flavor /, Q 2 is the invariant mass of a lepton pair, q^ (x, Q 2 ) and 
q^ A \x,Q 2 ) are respectively the quark and anti-quark distribution function with Bjorken 
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and by Gavin and Milanat 24 l One is Axi = a < L > A / E h , where a denotes the incident 
quark energy loss per unit length in nuclear matter, < L >a= 3/4(1. 2A 1//3 )fm is the average 
path length of the incident quark in the nucleus A. Another one is Ax 1 — f3 < L >\ / E h . 
Obviously, the quark energy loss is quadratic with the path length. In what follows, the two 
different parametrizations are called the linear and quadratic quark energy loss, respectively. 
The third form is Axi = kxiA 1 ^ 3 , which is named the incident-parton momentum fraction 
quark energy loss. In these three expressions, a, (3 and k can be extracted by a global analysis 
to nuclear Drell-Yan experimental data on the differential cross section ratio, respectively. 

The quark energy loss in target nucleus shifts the incident quark momentum fraction 
from x[ — x± + Ax i to x\ at the point of fusion. With adding the quark energy loss in the 
nucleus, the nuclear Drell-Yan differential cross section can be expressed as 

The parton distribution functions inside a nucleus from the Drell-Yan differential cross 
section have been found to differ notably from the corresponding ones in the free nucleon 
with the discovery of the nuclear EMC effect some twenty years ago(see Ref.[25], and ref- 
erences therein). Despite a significant worldwide effort in experiment and theory, there 
is as yet no consensus concerning the origin of this effect. In view of the importance for 
finding any new physical phenomena in the high-energy nuclear reactions, the global anal- 
yses of nuclear parton distribution functions, which parallel those for the free proton, have 
been performed in the past decade by different groups: HKM/HKN07[ 26 ' 27 ^, nDS' 28 ', and 
EPS09[ 29 l The four sets of nuclear parton distribution functions employed the existing ex- 
perimental data on nuclear structure functions from the electron and muon deep inelastic 
scattering. Unfortunately, the nuclear structure functions are composed of nuclear sea and 
valence quark distributions. The fact results in that the nuclear valence quark distributions 
are relatively well determined except for the small Bjorken variable x region, and nuclear 
antiquark distributions in small x region. It is difficult to constraint the antiquark distribu- 
tions at medium and large x region. It is expected that the nuclear Drell-Yan experimental 
data can pin down the nuclear valence quark distributions in the small x region and nuclear 
antiquark distributions in the medium x region 0.01 < x < 0.3. For this reason, HKN07 
and EPS09 added Fermilab E772 and E866 nuclear Drell-Yan data, and nDS included E772 
experimental data with difference from HKM. 
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FIG. 1: The nuclear modification of sea quark distribution at Q 2 = 50GeV 2 as a function of Bjorken 
variable x for tungsten nucleus. The solid, dashed, dotted and dash dot lines correspond to the 
the results given by HKM, HKN07, nDS and EPS09 nuclear parton distributions, respectively. 

The quantitative comparison between the different sets of nuclear parton distribution 
functions shows that the nuclear modification for valence quarks agrees nicely in the large- a; 
region x > 0.3. In other x region, HKN07, nDS and EPS09 give nuclear modifications 
relatively close to each other. Only the HKM displays a smaller antishadowing in the region 
0.01 < x < 0.3 than other sets, and no nuclear correction in x < 0.01. The nuclear 
modification of sea quark distribution for tungsten nucleus, RY(x,Q 2 = 50GeV 2 ) is shown 
in Fig.l, as a function of Bjorken variable x from the leading order HKM (solid line), HKN07 
(dashed line), nDS (dotted line), and EPS09 (dash dot line) nuclear effects. It is found that 
the nuclear modifications from different sets are relatively close to each other in the region 
0.01 < x < 0.08, however gives a clear deviation in the region x > 0.08. It is apparent 
that the nuclear modifications with including nuclear Drell-Yan data do not give a good 
consistency in the medium x region from HKN07, nDS and EPS09 parameterizations. The 
fact goes contrary to one's wishes. 
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3 The experimental data 

In our present analysis, the experimental data, providing the input for the value of the pa- 
rameter in the three representative expressions of quark energy loss, are taken from NA3^ 14 1 
and NA10 [151 Collaborations at CERN, and E772t 10 ] and E8661 17 ! Collaborations at Fermi- 
lab. The experimental data sets used here are summarised in Table I, in which the beam 
energy Eb eam of incident hadron, the projectile/target species, the covered domain on the 
momentum fraction of the hadron and target parton, and the number N of points in each 
data sample are specified. In total, our analysis has 269 data points, and 6 nuclei from 
Beryllium up to Platinum. 



TABLE I: Experimental data sets selected for the present analysis. 



Exp. E beam (GeV) Proj. 


Target 


x 1 


No. data 


X2 


No. data 


NAlOat 15 ] 140 


7T~ 


D, W 


0.39-0.82 


5 


0.163-0.360 


4 


NA3M 150 


7T~ 


H, Pt 


0.25-0.95 


8 


0.074-0.366 


7 


NAlObt 15 ] 286 


7T~ 


D, W 


0.22-0.83 


9 


0.125-0.451 


6 


E772I 16 ! 800 


P 


D, C, Ca, Fe, W 0.15-0.85 


122 


0.04-0.271 


36 


E866t 17 ] 800 


V 


Be, Fe, W 


0.21-0.95 


56 


0.01-0.12 


16 



To be emphasized, NA3 and NA10 data used cover the momentum fraction of the target 
parton from 0.074 to 0.366, and from 0.125 to 0.451, respectively. In the region x 2 > 0.1, the 
nuclear modification to sea quark distribution given by the different sets displays a gradually 
large deviation from each other with the increase of the momentum fraction of the target 
parton. Therefore, NA3 and NA10 data can show better the difference between HKM, 
HKN07, nDS and EPS09 parameterizations. Meanwhile, the bigger range coverage on beam 
energy and momentum fraction of the target parton can make us to find very clearly the 
quark energy loss effect on the nuclear Drell-Yan differential cross section ratio. 

4 Results and discussion 

In order to study the quark energy loss effect in the hadron-induced Drell-Yan reaction on 
nuclei, determine the values of the parameters a, (3 and k in quark energy loss expressions, 
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TABLE II: The x 2 /N-values computed using HKM, HKN07, nDS and EPS09 nuclear parton 
distribution functions without quark energy loss effect. The notation x\ and x 2 indicate the 
momentum fraction of the incident hadron and target parton,respectively. 



£jxp. uata 


tltS.Nl 


xlrvlN U 1 


nuo 




IN AlUa^xi ) 


no £!Ci 

zo.by 


QA 

0.04 


o.oi 


K HQ 

O.Uo 


IN AlUa(^2 J 


OK RQ 


K OR 


1U.4Z 


/i no 
4.yz 


AT A Q ( ™ \ 

IN K6\X\ ) 


k no 

o.yy 


A CI 

4.01 


4.40 


O.DO 


NA3(x 2 ) 


9.93 


6.98 


7.67 


6.35 


NAlOb(xi) 


3.81 


1.45 


1.60 


1.42 


NA10b(x 2 ) 


5.11 


0.65 


1.90 


0.86 


E772(xi) 


1.92 


1.41 


1.47 


1.33 


E772(x 2 ) 


4.83 


1.58 


1.79 


0.82 


E866(xi) 


1.44 


0.90 


1.19 


0.84 


E866(x 2 ) 


3.00 


1.17 


2.17 


0.95 



and investigate the dependence of quark energy loss on the nuclear parton distribution 
functions, we calculate in perturbative QCD leading order(LO) the Drell-Yan cross section 
ratio R*^J A2 on two different nuclear targets bombarded by hadron 

The comparison is performed with selected experimental data on the Drell-Yan differential 
cross section ratio. The integral range in above equation is obtained by means of the relative 
experimental kinematic region with neglecting the nuclear modifications in Deuterium. In 
our calculation as following, we use the four sets of leading order nuclear parton distribution 
functions together with CTEQ6L parton density in the proton I 30 ' , and parton density in the 
negative pion^ 31 !. 

When the incoming quark energy loss effect is neglected in the hadron-induced Drell- 
Yan reaction on nuclei, the calculated results are compared with the experimental data 
selected for our analysis on the Drell-Yan differential cross section ratio as a function of 
the momentum fraction of the incident hadron and target parton. The x 2 /N (N being the 
number of data points) computed are summarized in Table II by means of HKM, HKN07, 
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TABLE III: The values of a, (3 and k, x 2 /ndf and S factors extracted from each data sample with 
HKM nuclear corrections. The bottom row corresponds to the global fit of all selected data. 



Exp. data 


a ( X 2 /ndf, S) 


6 ( Y 2 /ndf, S) 




« ( X 2 /ndf, S) 




NAlOa(xi) 


1.53 ±0.14(1.61,1. 


.27) 


0.248 ± 0.022(1.49,1 


.22) 


0.0152 ± 


0.0015(4.12,2 


.03) 


NA10a(:c 2 ) 


1.46 ± 0.16(0.81,1. 


.00) 


0.230 ± 0.025(0.74,1 


.00) 


0.0210 ± 


0.0030(2.69,1 


64) 


NA3(xi) 


1.82 ±0.29(1.69,1. 


.30) 


0.310 ± 0.050(1.44,1 


.20) 


0.0120 ± 


0.0040(3.40,1 


84) 


NA3fx9) 


1.85 ± 0.36(2.53,1. 


.59) 


300 ± 060f2 24 1 


.50) 


0.0130 ± 


0.0050(5.78,2 


40) 


NAlObf^i) 


n ae 4- n 9nn ^1 1 

u.oo in u.zui _i_.o±,±. 


.14; 


140 ± 030C1 31 1 


14) 


n nn7n 4- 


n nm A7 i 


21) 


NA10b(z 2 ) 


0.79 ±0.14(0.45,1. 


.00) 


0.125 ±0.023(0.45,1 


.00) 


0.0083 ± 


0.0016(0.88,1 


.00) 


E772(xi) 


1.26 ±0.16(1.33,1. 


,15) 


0.230 ±0.030(1.35,1 


.16) 


0.0042 ± 


0.0005(1.26,1 


.12) 


E772(x 2 ) 


1.29 ±0.13(1.10,1. 


.05) 


0.230 ±0.030(1.53,1 


.24) 


0.0066 ± 


0.0006(1.69,1 


.30) 


E866(xi) 


1.28 ±0.22(0.79,1. 


.00) 


0.190 ±0.030(0.80,1 


.00) 


0.0026 ± 


0.0004(0.77,1 


.00) 


E866(x 2 ) 


1.27 ±0.23(1.12,1. 


.06) 


0.190 ±0.030(1.13,1 


.06) 


0.0035 ± 


0.0006(0.78,1 


.00) 


Global fit 


1.21 ±0.09(1.07,1. 


.03) 


0.190 ±0.020(1.11,1 


.05) 


0.0037 ± 


0.0003(1.21,1 


.10) 



nDS and EPS09 nuclear parton distribution functions, respectively. Usually, if the x 2 /N 
is not much larger than one, the theoretical results are considered as being statistically 
consistent with the experimental data. As seen from the Table II, each analysis on nuclear 
effects has consistently much larger x 2 /N values for NA3 and NAlOa data from the negative 
pion incident Drell-Yan reaction on nuclei. Interestingly, very larger x 2 /N value is 23.69 and 
25.68 for NAlOa(xi) and NAl0a(x 2 ) data, respectively, from the HKM nuclear effects. 

To determinate the optimal parameter from each experimental data set, we adopt the % 2 
analysis method described in Ref.[10]. The obtained results are summarized in Table III. by 
combining the HKM cubic type of nuclear parton distributions with linear, quadratic and 
incident-parton momentum fraction quark energy loss, respectively. The values of a, (3, k 
extracted from the individual fit of each data sample, as well as their corresponding rescaled 
error, x 2 P er number of degrees of freedom (x 2 /ndf) and S factors, are listed in Table III, in 
which the bottom row corresponds to the global fit of all selected data. As can be found from 
Table III compared with Table II, x 2 values with quark energy loss effect, overall, are much 
smaller than those with only HKM nuclear effects of parton distributions. The agreement 
of theoretical calculations with experimental data has a significant improvement. However, 
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FIG. 2: The nuclear Drell-Yan cross section ratios Ra 1 /a 2 ( x i) an d Ra 1 /a 2 (. x 2) by using HKM 
nuclear parton distributions. The solid curves correspond to the results with only nuclear effects 
of parton distributions. The dashed, dotted and dash dot curves show the combination of nuclear 
effects of parton distributions with linear, quadratic and incident-parton momentum fraction quark 
energy loss, respectively. The relative optimal parameter is taken from the fit to corresponding 
data sample. The experimental data are taken from NA3 [14] and NAl0a [151 . 

the computed results with the incident-parton momentum fraction quark energy loss have 
yet a significant deviation from NA3 and NAlOa data sets. 

Regarding the linear, quadratic and the incident-parton momentum fraction quark energy 
loss, the global fit of all data makes a = 1.21 ± 0.09 with the relative uncertainty 5a/ a ~ 
7% and x 2 /ndf = 1.07, f3 = 0.19 ± 0.02 with 5/3//3 ~ 10% and x 2 /ndf = 1.11, and 
K = 0.0037 ± 0.0003 with Sk/k ~ 8% and x 2 '/ndf = 1.21. 

To demonstrate intuitively the energy loss effect of an incoming quark on the nuclear 
Drell-Yan cross section ratio, the calculated results combining HKM cubic type of nuclear 
parton distributions are compared with NAlOa and NA3 data in Fig. 2. It is necessary to 
note that NA3 Collaboration provided the negative pion incident Drell-Yan cross section 
ratio Ru/pt as a function of parton momentum fraction. If the ratio R,H/pt is transformed 
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TABLE IV: The values of a, (3 and k with % 2 /ndf and S factors extracted from the global fit of 
all data by using HKN07, nDS and EPS09 nuclear corrections, respectively. 





a HKN07 ( x 2/ ndf: S ) 


a nDS ( x 2 /ndL g) 


a SP509 ( x 2/ ndf: S ) 


Global fit 


0.64 ±0.09(1.09,1.04) 


0.73 ±0.09(1.08,1.04) 


0.23 ±0.07(1.05,1.02) 




P HKN07 ( x 2 /ndL S) 


pnDS ( x 2 jndL S) 


/?™ ( x^/ndf, S) 


Global fit 


0.103 ±0.014(1.09,1.04) 


0.122 ±0.013(1.08,1.04) 


0.042 ±0.015(1.05,1.02) 




K HKN07 ( x 2/ ndfj S ) 


K nDS ( x 2 /ndfi g) 


K EP509 ( x 2 /ndL S) 


Global fit 


0.0023 ±0.0003(1.08,1.04) 


0.0026 ±0.0004(1.08,1.04) 


0.0009 ± 0.0004(1.05,1.02) 



to Rpt/n, the ratio Rpt/H has the same tendency as Rw/d at the beam energy 140 GeV. 
It can be seen from the Fig. 2 that the calculated differential cross section ratios are nearly 
same from the linear and quadratic quark energy loss. The theoretical prediction from the 
incident-parton momentum fraction quark energy loss is not in agreement with the NA3, 
NAl0a(x2) and NAl0a(x 1 ) data in x\ < 0.4. Therefore, we can conclude apparently that 
the existing experimental data from lower incident beam energy rule out the possibility of 
the incident-parton momentum fraction quark energy loss. Whether the quark energy loss 
is linear or quadratic with the path length is not determined. 

In order to quantify the sensitivity of our computed results on quark energy loss with 
respect to the strength of nuclear corrections to the parton distributions, the global fit 
analysis is carried out by using HKN07, nDS and EPS09 parameterizations. The results 
given from the three parametrizations of quark energy loss are summarized in Table IV with 
the values of a, (3 and k, their corresponding rescaled error, x 2 jndf and S factors. It can 
be found that the obtained parameter values in three expressions for quark energy loss are 
smaller than that by HKM nuclear effects. It directly reflects the deviation between HKM 
nuclear corrects to sea quark distribution and other three sets. The calculated results are 
compared with NA3' 14 ' and NAlOa' 15 ' data including quark energy loss effect and nuclear 
effects of parton distributions from HKN07, nDS and EPS09 sets in Fig. 3, Fig.4 and Fig. 5, 
respectively. It is shown from these figures that the computed results from the linear quark 
energy loss are identical to those from the quadratic quark energy loss. The theoretical 
prediction on Rw/d from the incident-parton momentum fraction quark energy loss exists a 
significantly large deviation from that by the linear (or quadratic) quark energy loss in the 
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X 2 X-i 

FIG. 3: The nuclear Drell-Yan cross section ratios Ra 1 /a 2 ( x i) an d Ra 1 /a 2 { x i) by using HKN07 
nuclear parton distributions. The other comments are the same as those in Fig. 2. 

region x 2 > 0.3 and X\ < 0.4. The tendency of Rw/d as a function of quark momentum 
fraction does not support the possibility of the incident-parton momentum fraction quark 
energy loss. Additionally, it can be found that with independence on the nuclear modification 
of parton distributions, NA3 experiment rules apparently out the incident-parton momentum 
fraction quark energy loss. In view of the large experimental error in NA3 and NA10 data, it 
is desirable to operate precise measurements on the nuclear Drell-Yan reactions from lower 
incident beam energy. 

5 Summary and concluding remarks 

In summary, the available data on nuclear Drell-Yan differential cross section ratio as a 
function of the quark momentum fraction have been analysed with three parametrizations of 
quark energy loss and four typical sets of nuclear parton distribution functions. It is found 
that with independence on the nuclear modifications of parton distributions, the experi- 
mental data from lower incident beam energy rule out the possibility of the incident-parton 
momentum fraction quark energy loss. The existing experimental data do not distinguish 
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FIG. 4: The nuclear Drell-Yan cross section ratios Ra 1 /a 2 ( x i) an d Ra 1 /a 2 ( x 2) by using nDS 
nuclear parton distributions. The other comments are the same as those in Fig. 2. 

between the linear and quadratic dependence of quark energy loss. It is worth to mention 
that the mean energy loss is employed in our calculations. In hot and dense matter, however, 
the mean energy loss of the highly energetic partons would be considered very simplistic. 
Rather, it is now accepted that at least the probability distribution P(AE, L) of energy loss 
AE given a path L is the relevant quantity, which then needs to be averaged over geometry, 
i.e. a calculation needs to include explicitly both dynamical fluctuations given the same 
path, and fluctuations of the path a quark takes through the medium. It is possible that the 
distinction between linear and quadratic energy loss as observed in present analysis is lost 
by not accounting for fluctuations. However, while the quadratic dependence of energy loss 
is argued to arise from the Landau- Pomeranchuk-Migdal effect, it is now known that effec- 
tively due to finite energy corrections even an Landau-Pomeranchuk-Migdal-driven radiative 
energy loss reverts to an approximately linear dependence quickly' 32 ' 33 ' . 

From the global fit of all selected data, we obtain the quark energy loss per unit path 
length a = 1.21±0.09 GeV/fm by HKM nuclear parton distribution functions. By combining 
our previous discussion on the semi-inclusive deep inelastic scattering of lepton on nuclear 
targets' 10 ', our result on the mean energy loss per unit length of an incoming quark is 
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FIG. 5: The nuclear Drell-Yan cross section ratios Ra 1 /a 2 ( x i) an d Ra 1 /a 2 { x 2) by using EPS09 
nuclear parton distributions. The other comments are the same as those in Fig. 2. 

not in support of the theoretical prediction: the mean energy loss of an outgoing quark is 
three times larger than that of an incoming quark approaching the medium ' 34 l. In ultra- 
relativistic heavy-ion collisions, however, the medium-modified fragmentation function is 
obtained from a computation of an in-medium parton-shower followed by hadronization' 35 ' 36 ', 
i.e. what matters for the final state is not the energy loss of a single quark but rather 
the modified development of a parton shower. This virtuality evolution may explain the 
difference between the energy loss of incoming and outgoing quarks - while incoming quarks 
all in all probably can be considered on-shell, outgoing quarks are significantly off-shell due 
to the hard scattering. 

In addition, the obtained value of the parameter in the quark energy loss expression from 
HKM nuclear effects is larger than that using other three sets of nuclear parton distribution. 
Our calaulated results show that the energy loss effect of the incident quark has a distinct 
impact on the Drell-Yan cross section. It directly brings about an overestimation for nuclear 
correct to the sea quark distribution if leaving the quark energy loss effect out. Besides, 
the Drell-Yan single differential cross section as a function of the target parton momentum 
fraction is dominated by nuclear sea and valance quark distribution, which is similar to the 
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nuclear structure function in charged-lepton deep inelastic scattering on nuclei. In order to 
make the flavor decomposition of nuclear parton distribution functions, we need to resort to 
the neutrino deep inelastic scattering data. Several works have studied the nuclear effects 
in the neutrino-nucleus charged- current inelastic scattering process' 37 ' 38 ' 39 '. We suggest that 
the new global analysis of nuclear parton distribution functions should employ the available 
experimental data on structure function from neutrino and charged-lepton deep inelastic 
scattering on nuclei. 

It is worth noting that the fractional energy loss does not provide a good description of the 
data in hot-dense matter physics' 40 ] because the fast parton propagation in cold nuclear and 
hot-dense matter contains different physics' 41,42 ]. Moreover, in our present work, the used 
experimental data from NA10 Collaboration recorded the x 2 dependence of nuclear Drell-Yan 
cross section ratio from 0.12 < x 2 < 0.45, which can be covered by Fermilab E906/SeaQuest 
experiment' 43 ^. Therefore, we desire that our paper can provide useful reference for E906's 
insight on the enengy loss of an incoming quark propagating in cold nucleus. 
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